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ABSTRACT 

We study the effect of primordial isocurvature perturbations on non-Gaussian proper- 
ties of CMB temperature anisotropics. We consider generic forms of the non-linearity of 
isocurvature perturbations which can be applied to a wide range of theoretical models. 
We derive analytical expressions for the bispectrum and the Minkowski Functionals 
for CMB temperature fluctuations to describe the non-Gaussianity from isocurva- 
ture perturbations. We find that the isocurvature non-Gaussianity in the quadratic 
isocurvature model, where the isocurvature perturbation S is written as a quadratic 
function of the Gaussian variable a, S = a 2 — (a 2 ), can give the same signal-to-noise as 
/nl = 30 even if we impose the current observational limit on the fraction of isocurva- 
ture perturbations contained in the primordial power spectrum a. We give constraints 
on isocurvature non-Gaussianity from Minkowski Functionals using the WMAP 5- 
year data. We do not find a significant signal of isocurvature non-Gaussianity. For the 
quadratic isocurvature model, we obtain a stringent upper limit on the isocurvature 
fraction a < 0.070 (95% CL) for a scale invariant spectrum which is comparable to 
the limit obtained from the power spectrum. 

Key words: Cosmology: early Universe - cosmic microwave background - methods: 
statistical - analytical 



1 INTRODUCTION 

Recent observational progress in cosmology represented 
by surveys such as WMAP and SDSS has enabled a 
detailed analysis of cosmic density fields to investigate 
the physics of the early universe. In particular, (non- 
)Gaussianity of primordial density fields has been received 
much attention recently as a key observational probe to 
differentiate models of the early universe. In the simplest 
single field inflationary scenario, quantum fluctuations 
of the inflaton during inflation are assumed to be the 
origin of cosmic density fluctuations and such a model 
predicts adiabatic and almost Gaussian primordial fluctu- 
ations. However, other generation mechanisms of density 
fluctuations such as the curvaton scenario ( Mollcrach 
Il99d : iLinde fc Mukhanovl 1 19971 : iMoroi fc Takahashil l200lT 
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lEnqvist fc Slothl 120021: iLyth fc Wandsl l2002lh modulated 



rehea ting (|Kofman 20031 : Dvali. Gruzinov fc Zaldarriagal 



l2004h and so on have been proposed. In these mechanism, 
the nature of primordial density fluctuations can be very 
different from that of a simple inflation model, in partic- 
ular, in terms of Gaussianity of the fluctuations. In fact, 
it has been shown that large no n-Gaussianity can be gen- 
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iLehners fc Steinhardtl |2008|). single field inflation wit h a 
feature in its potential ( Chen. Easther fc Liml l2007al ) . a 
Gaussian-sq uared component in curva t ure or entropy per- 
turbations jLinde fc Mukhanoy 1 19971 ; iBoubekeur fc Lvthl 
120061; Suvama fc Takahashi 20081) a nd multi-brid inflation 
(|Sasakill2008l : iNaruko fc Sasakill2009l ). 

The primordial non-Gaussianity has been quantita- 
tively measured by higher-order statistics such as bispectra 
and Minkowski Func tionals from CMB temperature maps 
obtained by WMAP (IKomatsu et al.ll200Sl: ICreminelli et all 



120071 ; lYadav fc Wandeltll2008l; iHikage et alJbOOsh and also 
from large-scale structure (jSlosar et al.ll2006 ) . So far the ob- 
servational results are consistent with the Gaussian hypoth- 
esis. There is, ho wever, a hint of primordial non-Gaussianity 
at 2-3 a level l|Yadav fc WandehJ 120081 : iKomatsu et al] 
120081 ). Since the non-Gaussianity predicted in the simplest 
inflation model is too small to be detected by current obser- 
vations, if a non-Gaussian signal is observed and it originates 
from primordial fluctuations, other mechanisms beyond the 
simplest model would be required in the dynamics of early 
universe. 

As another probe of the early universe, the adiabatic- 
ity of primordial density fields has also been the sub- 
ject of intense study. In fact, current cosmological obser- 
vations of TT and TE spectra of CMB with some other 
distance measurements such as type la supernovae (SNe) 
and baryon acoustic oscillation (BAO) have already con- 
strained the fracti on of isocurvature perturbations to be 
less than 10% (e.g., Bean et al . 2006; K awasaki fc Sekiguchil 
120071 ; IKomatsu et al.ll2008l ). Examples of isocurvature per- 
turbations along with non-Gaussianity have been dis- 
cussed in the context of non- Gaussian field potentials 
dLinde fc Mukhanovll 19971; iPeeblesil 19991 ; IBoubekeur fc Lvthl 
120061 ; ISuvama fc Takahashil |200Sh. the curvaton sce- 
nario dLyth. Ungarelli fc Wandsl 120031; iBartolo et al] l20oi 
iBeltranl 120081; iMoroi fc Takahashil I200ST ), modulated re- 
hea ting (IBoubekeur fc Creminelh"l2006l), ba ryon asymme- 
trv dKawasaki. Nakavama fc Takahashill2009l ) and the axion 
l|Kawasaki et al]|2008l . 120091 '). In particular, non-Gaussianity 
generated from isocurvature pertur bations has been sys- 



tematically investigated recently dKawasaki et ah I l200Sl: 



lLanglois. Vernizzi fc Wandsll2008b1 ; iKawasaki et al]|2009f ) 

In this paper we discuss non-Gaussianity generated 
from the non-linearity of isocurvature perturbations and 
study a constraint on non-Gaussianity from isocurvature 
perturbations using Minkowski Functionals. For this pur- 
pose, we derive theoretical expressions for bispectra and 
Minkowski Functionals to characterize non-Gaussianity in 
CMB temperature maps generated from primordial mixed 
perturbations of adiabatic and isocurvature components. We 
characterize the non-linearity of isocurvature perturbations 
in two different forms which are theoretically motivated; one 
is a Gaussian variable plus its quadratic correction (Linear 
Model) . The other form is given as a quadratic of Gaussian 
variables without a liner term (Quadratic Model). These two 
generic forms are applicable to a wide range of isocurva- 
ture models listed above. Then we give actual limits on the 
isocurvature non-Gaussianity using the WMAP 5-year data. 
As far as we know, this is the first attempt to put a limit 
on isocurvature perturbations from the non-Gaussianity of 
CMB anisotropies. 

In this paper, we focus on CDM isocurvature perturba- 



tions. However, it is straightforward to apply our method to 
other types of isocurvature perturbations including baryon 
and neutrino ones. We adopt a set of cosmological pa- 
rameters at the maximum likelihood values for a power- 
law AC DM model obtained from the WMAP 5- year data 
only fit (|Dunklev et alj|2009h ; Sl b = 0.0432, fi cdm = 0.206, 
Oa = 0.7508, H = 72.4 kms" 1 Mpc" 1 , r = 0.089, and 
?!</, = 0.961. The total amplitude of primordial power spec- 
tra is set to be A to t(fc = 0.002Mpc _1 ) = 2.41 x 10~ 9 . 

This paper is organized as follows; in &\ we give two 
different forms of the non-linear isocurvature perturbations 
called the "Linear Model" and the "Quadratic Model". In 
we derive analytical expressions for bispectra to de- 
scribe isocurvature non-Gaussianity in CMB temperature 
anisotropies in the Linear Model. The isocurvature non- 
Gaussianity in the Quadratic Model is described in 331 In 
[J5] we present generic perturbative formulae of Minkowski 
Functionals that can be applied to CMB temperature maps 
with adiabatic and isocurvature non-Gaussianity. In 
we give limits on isocurvature non-Gaussianity from the 
WMAP 5-year temperature maps using Minkowski Func- 
tionals. In 53 we discuss implications of our results for an 
axion isocurvature model. |8]is devoted to a summary and 
our conclusions. 



2 NON-LINEAR ADIABATIC AND 
ISOCURVATURE PERTURBATIONS 

We consider the admixture of an adiabatic perturbation £ 
with a CDM isocurvature perturbation 5. The curvature 
perturbation is written up to second order in a local form 

a.s 



C = 0+ jJkW -fa », 



(1) 



where <fi is the linear term of ( that obeys a Gaussian statis- 
tics. The non-linear parameter /nl represents the quadratic 
amplitu de of the curvature pertur bation <E> during the mat- 
ter era (|Komatsu fc Sper gcl 2001), which is related to £ by 
$ = (3/5)C 

An isocurvature perturbation 5 between matter and ra- 
diation is defined as 



5 = 



5p m 3Jp 7 



4p 7 



(2) 



where p m is the matter energy density and p 7 is radiation 
energy density. In this paper we focus on a CDM isocurva- 
ture perturbation. 

We consider the non-linearity of the isocurvature per- 
turbation in two different forms. One is a local form similar 
to the equation ([T]) which has a linear (Gaussian) term with 
a quadratic correction 



I. Linear Model 



(n)), 



(3) 



where rj is a Gaussian variable and its non-linearity is 
characterized by /^°'. The isocurvature non-Gaussianity 
in t his form was s t udied in the context of the ax- 
ion (Kawasaki ct al. I I2008D and the curvaton scenario 
(|Langlois. Vernizzi fc Wandj|2008bh . 

The other is the case where the linear term 
is negligible compared with the quadratic term (e.g., 
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Linde fc Mukhanovlll997l; |Peebleslll999l ; IBoubekeur fc Lvthl 
20061 ; iKawasaki et al.ll2008h ; 



A^k). 



(12) 



(4) 



II. Quadratic Model : 5 = 

where a obeys Gaussian statistics. iLinde fc Mukhanovl 
\ 19971 ) proposed a scenario to generate the quadratic form 
of isocurvature perturbations by introducing a massive free 
scalar field oscillating around the vacuum state which has 
a zero value. In this scenario, the isocurvature fluctuation 
has a blue spectrum and thus we here consider a wide range 
for the spectral index rang i ng fro m 1 to 3 in the Quadratic 
Model. IBoubekeur fc Lvthl (|2006l ) showed that a Gaussian- 
squared component of the primordial curvature perturbation 
would be bounded at 10% level by the WMAP bound on the 
bispectrum. 

The auto and cross-correlation power spectra for fluc- 
tuations of / and J are defined as 



(J k J k ,) = (27r) 3 <5< 3 3) (k + k')P7./W, 

and then its dimensionless power is given by 

Pjj(k)k 3 
Alj{k) = 2^ ' 
where / and J denote ( or 5. 

We assume the following power-law form of auto and 
cross power spectra for the Gaussian variables <j> (eq. [I]), 77 
(eq-GS), an d o (eq.gj); 



(5) 



(6) 



k \ nx — 1 



A xx (k) = 
A XY (k) = (Ax A Y f 1 2 cos 6 XY (A) 



{n x +n Y )/2-l 



(7) 



(8) 



where X and Y denote cj>, 77 or a. As /nlA^ is observa- 
tionally limited to be much smaller than unity, the power 
spectrum of the primordial adiabatic perturbation is given 
by 



A fC (fc)~ A^(fe). 



(9) 



We define the fraction of isocurvature perturbation as 
Pss(ko) 

~ P a (k ) + Pss(k )' { ' 

which is the same definition as Bean et al.1 <|2006h for ex- 
ample. We here set ko = 0.002Mpc _1 . The parameter a 
is related to another c ommon parameter of adiabaticity 
S id^ 0-S-> eq- ( 41 ) of iKomatsu et all I2008T ) as = 
[a/(l - a)] 1/2 /3- The upper limits of a from WMAP, BAO 
and SN combined are given by 0.067 (95% CL) for axion- 
type (cosO^s — 0) and 0.0037 (95% CL) fo r curvaton-type 
(cos P es = — 1) isocurvature perturbations (|Komatsu et al.l 
l2008h . 



3 NON-GAUSSIANITY OF ISOCURVATURE 
PERTURBATIONS I: LINEAR MODEL 

3.1 Initial Perturbation 

In the Linear Model (eq.[3]), power spectra for the isocur- 
vature perturbation and its cross term with the adiabatic 
perturbation become 



Here, we have used the fact that < 1/y/Ass — 

aA^ ~ 10 6 . In the case that J^ ' is larger than this 

upper limit, the linear term of 77 in Eq. ((3j) is negligible so 

that the model should be described by the quadratic model. 
The ratio of the amplitude of the power spectra between 

<j> and 77 is written in terms of a (eq.[l0]) as 



1-a 



(13) 



We define the bispectra of £, S and their mixed contri- 
bution as 



</ kl Jk 2 ^k 3 ) = (27T) 3 5g ) (k 1 +k 2 +k,) 

xBij K {ki,k2, k 3 ), 



(14) 



A ss (k) 



Ar/r](k), 



(11) 



where I, J, and K denote £ or S. 

When 4> an d 77 are initially uncorrelated, the adiabatic 
and isocurvature modes evolve independently under the lin- 
ear approximation. There exist bispectra only from each 
mode given as 

Bac(ki, k 2 , k 3 ) = ^fNL[P4,4,(k 1 )P 4 , <l ,(k2) 

+P^(k 2 )P^(k 3 ) + P^(k 3 )P^(ki)], 

(15) 

Bsss(ki,k 2 ,k 3 ) = 2/^° ) [P^(fci)P m (fc>) 

+P vv (k 2 )P vv (k :i ) + P w (fe 3 )P w (fti)]- 

(16) 

When 4> and 77 are initially correlated, the following cross- 
correlation terms may become important 

B( (( s)(ki,k 2 ,k 3 ) = B QC s{ki,k 2 ,k 3 ) + B Q sc{k 1 ,k 2 ,k 3 ) 
+B S (;c(ki,k 2 , k 3 ) 
= 2f^° ) [P^(k 1 )P^(k 2 ) 

+P^(fo)P^(fe) + P^v(k'i)P^(ki)} 
+ ^/nl[P^(A;i) {P^(fo) + P^(fcs)} 

+P^(fo) {P^(fcs) + P^(fci)} 
+P^(fc 3 ){P^(fci) +P^(fe)}], (17) 
B( ( ss)(ki,k 2 ,k 3 ) = P ? ss(fci, fc 2 , fc 3 ) + P sc<s (fci, fc 2 , fc 3 ) 
+B S sc(ki,k 2 ,k 3 ) 

= ^/NL[P^(fcl)P^(fe) 

+P^{k 2 )P <t>v {k 3 ) + P^(fc 3 )P^ ) (fci)] 
+2f^° ) [P^(k 1 ) {P vv (k 2 ) + P vv (k 3 )} 
+P^(fe) {P w (fci) + Pvv(ki)} 
+P^{k 3 ) {P vv {ki) + P vv (k 2 )}]. (18) 



3.2 Bispectra of CMB Temperature Anisotropy 

Harmonic coefficients a; m of CMB temperature anisotropics 
AT IT are defined as 

(19) 



ai m = I dn-^(n)Y^(n 
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Figure 1. Angular power spectra C; for adiabatic and isocur- 
vature perturbations SS. The plotted adiabatic perturbation has 
the spectral index = 0.96 (solid). For isocurvature perturba- 
tions, the spectral index of a Gaussian variable is 1 (long- 
dashed) in the Linear Model and n a is 1 (short-dashed), 1.5 
(dotted) and 2 (dot-dashed) in the Quadratic Model (see jj4}. 
The isocurvature fraction a is set to be 0.067 (n v = 1 and 
n a = 1), 0.008 (n CT = 1.5), and 0.001 (n CT = 2) defined at 
fc = 0.002Mpc _1 . 



Introducing the radiative transfer function g^ { (or g^i), they 
are related to £ (or S) as 



a l n 



s 

(II, 



I ° 
T a lr, 



4tt(- 



= 4 7 r(-i) i 



ri 3 k 

(2tt) 3 

d 3 k 



Ck^(*)yim(i), 



? 5 kff ? i (fc)y i * m (k). 



(20) 
(21) 

(22) 



(2tt) 3 

The angular power spectra of adiabatic and isocurvature 
components are 



c 



= - / k 2 dkP <c (k)g ( Tl (k) 2 , 



k 2 dkP SS (k)gUkY 



(23) 



(24) 



Fig. Q] shows the angular power spectrum from adia- 
batic and isocurvature perturbations. The radiation trans- 
fer functions for adiabatic and isocurvature perturbations 
are computed using the pub licly-available CMBFAST code 
l|Seliak fc Zafdarriagal Il996t ). The spectral index of isocur- 
vature perturbations in the Linear Model is n n = 1. We set 
a = 0.067 for n v = 1 at k = 0.002Mpc~ 1 which is a 2a 
uppe r limit from the WMAP 5-year paper (jKomatsu et alj 
l200Sh . 

The total angular bispectrum of CMB is written as the 
sum of bispectra with different combinations of adiabatic 
and isocurvature components: 



— b hhh> 

UK 



where each component is defined as 

)/ — y hl 2 l~ " 



, I .1 K 

1 m l 

a l 2 m 2 a l 371 



mi m,2 rlJK 



(25) 



(26) 



dh Y hmi (n)Yi 2m2 (h)Yi 37n3 (h) 



(2/i + l)(2Z 2 + l)(2Z 3 + l) 



4tt 







h 




h 



h 
m 2 



h 
m 3 



(28) 



where I and J denote £ or S. Then, 6{^j 3 j 3 can be written 



kfdki / k\dk 2 / k 2 dk 3 tor (&i).7ii (for) 



,ijk 8 / 2 j 
b hhh = — I r dr 



XffTi 2 (fa) ji 2 (k2r)g T i 3 (k 3 )ji 3 {k3r)BijK(ki,k2,k 3 ). (29) 

The explicit form of adiabatic, isocurvature bispectra and 
their cross-correlations are given in Appendix lAl 

Fig. [2] shows each component of the CMB bispectrum 
for equilateral triangles (h = h = h ~ I)- The purely 
adiabatic component bj^ with /nl = 50 is plotted with 
a solid line. A long-dashed line represents a mixed compo- 
nent ftjjj with | cosf^l = 1 and the curvaton-type upper 

( CSS) 

limit a = 0.0037. The other mixed component 6^ with 
j cos#0 ?) ] =0.1 and the axion-type limit a — 0.067 is plotted 
with a short-dashed line. The pure isocurvature component 
b?if S with the axion-type limit a = 0.067 is plotted with 
dotted lines. The non-linearity of isocurvature perturbation 



.(ISO) 
/NL 

perturbation is 1 



/nl"^ ^ s se t t° be 10 4 . The spectral index of isocurvature 



We numerically estimate how large / 



(ISO) 



the non- 



Gaussianity from isocurvature perturbations, should be 
to obtain the same values of the signal-to-noise ratio of 
the CMB bispectrum as the one derived from the non- 
Gaussianity from purely adiabatic perturbations character- 
ized by /nl. The signal to noise ratio is denned as 



2<!i<! 2 <! 3 



r 2 



UK 



Ci 1 Ci 2 Ci 3 A 



hl 2 l3 



(30) 



where Ci is C; of purely adiabatic perturbations (a = 0) 
including noise and Ii x i 2 i 3 is defined as 



'hhh 



(2Zi + l)(2Z 2 + l)(2Z 3 + l) 



4tt 







h 





.(31) 



The factor Aj^jg is equal to 6 (h = l 2 = h), 2 (h — l 2 
or h = h), and 1 (h / l 2 / l 3 ). WMAP beam window 
functions are included both in the bispectrum and in the 
power spectra. The homogeneous noise distribution for the 
WMAP 5-year data is used to estimate the noise. 

The pure isocurvature term b sss with n v = 1 has non- 
Gaussianity corresponding to 



/NL " 15 (aS7) 



.(ISO) 
iNL 

10 4 



(32) 



where the above equations are valid when a 1. It is found 
that isocurvature non-Gaussianity reaches /nl ~ 10 only 
if the non-linearity in isocurvature perturbations /nl i s 
of order 10 4 because isocurvature non-Gaussianity in the 
Linear Model is suppressed by a 2 . 

When there is a weak correlation between <j> and r\ 
(|cos<V,| < 1), b< css > has non-Gaussianity corresponding 
to 



(27) /nl =± 9.7 



f— 

V 0.067 



3/2 



0.1 



.(ISO) 
■'NL 

10 4 



(33) 
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Figure 2. CMB angular bispectra of equilateral configurations 
l 2 (l + l) 2 6;n/(27r) 2 in the Linear Model; the adiabatic compo- 
nent (solid), the mixed components biff S ^ (long-dashed) and 
b^f 3 * 1 (short-dashed) and the isocurvature component bf^S (dot- 
ted). Upper (Lower) panel shows the positive (negative) side of 
bispectra plotted in logarithmic scale. The adiabatic perturba- 
tion has the power-law index = 0.96 and its quadratic ampli- 
tude /iml = 50. The isocurvature perturbation has n v = I and 



J i 



(ISO) 



10 4 . The fraction of the isocurvature power spectrum a 



is set to be an axion-type upper limit 0.067 for bf^ s and b^^ 3 " 1 
with a weak correlation | cos 8^ | =0.1 and a curvaton-type upper 



limit 0.0037 for b 



(CCS) 



with a strong correlation | cos i 



<,»>>/ 



If the correlation is strong | cos 9^ ~ 1| (e.g., curvaton-type 
isocurvature perturbation), the other correlated term b^^ s ' 
becomes important with a more severe limit a < 0.0037 for 
curvaton-type isocurvature perturbations; 



h 



7.2 



\ 0.0037/ 



'01/ 



i(ISO) 
J NL 

10 4 



(34) 



In both cases, however, isocurvature perturbations need to 
have a strong non-linearity f^°^ ~ 10 4 in order to gen- 
erate non-Gaussianity corresponding to /nl ~ 10. Our 
result is consistent with previous theoretical estimation s 
jKawasaki et al.ll200Sl ; lLanglois, Vernizzi fc Wandsll2008bh . 



where a finite box-size L gives an infrared cutoff. To avoid 
assumptions at scales far beyond the present h orizon H^ 1 , 
L should be set not too much bigger than H^ 1 (|Lvthll2007h . 
Hereafter we set L — 30Gpc. 

Using the power-law form as in the equations (0, the 
isocurvature power spectra are written as 



/ k \2(n<x-l) 

A S s(k) = AlF(n a ) (^J 

where the factor F is approximately given by 



(37) 



F(n s _ / 41og(fci) ( n<7 = 1), 

F ^)-\ A(n a {kL)^) K/l), (38) 

As we take L to be much larger than the range of scales we 
are interested in, the k dependence in F is weak. 

As 5 is the square of a Gaussian variable, the cross- 
correlation coefficient becomes nearly zero regardless of 
cos 6>0 CT ; 

(A CC A S5 )V2 



^ /nl(A^F) 1/2 COS 2 Ofra 

< /nl(^F) 1/2 <G>(TT 2 ). 



(39) 

Following the definition of a in the equation (|10[) . the 
ratio between A$ and A a becomes 



1/2 



-1/2 



(40) 



A^ " 

The amplitude of a has an additional factor (A^F)" 1 ^ 2 ~ 
O(10 4 ) relative to the amplitude of (f> because the isocurva- 
ture perturbation is given as quadratic in a (eq. U). 

The bispectra from isocurvature perturbations are writ- 
ten as 

P„(P) 



B SS s(ki,k 2 ,k 3 ) = - . . . 3 
Jl-i \ An ) 

xIP^dkr+pDP^dka-pl) 
+P CT a(|k 2 + p|)P CTCT (|k3-p|) 
+P CTCT (|k 3 + p|)P CTCT (|ki-p|)]. (41) 

When we extract the dominant contributions around the 
oles, the equation (|41|l is approximately written as 
Kawasaki et al.ll200&f ) 

B SS s(ki,k 2 ,h i ) ~ 2A tTa (k b )F(n rT )[P tTa (k 1 )P tT{T (k 2 ) 

+P«<r(k 2 )P rTtT (k 3 ) + P CTCT (fe)Prcr(fel)], 

(42) 

where kb =Min(fci, k 2 , k 3 ). 

The cross-correlation terms arc 



4 NON-GAUSSIANITY OF ISOCURVATURE 
PERTURBATIONS II: QUADRATIC MODEL 

4.1 Initial Perturbation 

In the Quadratic Model (eq.[4]), the auto and cross power 
spectra of £ and 5 become 



d 3 p 



(35) 



Pss(k) = 2 J j^P„(p)Pa„(\k + p\. 

|/nl I i ^_P 0CT (p)P 0CT (|k + P |), (36) 



P (ccs) (fci,fc 2 ,fc 3 ) = 2[P^{k 1 )P^{k 2 ) + 

P^(fc2)P^(fe 3 ) + P^(k 3 )P^{ki)], 

(43) 



24 



B( ( ss)(ki,k 2 ,k 3 ) = — /nl 



<f j p 

-1 



p CTCT (p) 



K[P^ ff (|ki-|-p|)P0 CT (|k 2 -p|) 
+P^(|k 2 -r-p|)P^(|k 3 -p|) 
+P^(|k 3 -r-p|)P^(|ki-p|)]. (44) 

The amplitude of each bispectrum component relative 
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to the pure adiabatic one at ko is estimated as 



Bsss 
B ccc 



B 



(CCS) 



B, 



ccc 



B 



(CSS) 



B 



5 
3 
5 
3 

= 3 



/nl 
/nl 



1 - a 



3/2 



(^) 



-1/2 



CCC 



1 



1/2 



70(J ■ 



-1/2 



(45) 
(46) 
(47) 



The bispectra B^qs) and Bsss has an additional factor 
(A^F)- 1 ^ ~ O(10 4 ) relative to B ccc . This is because the 
isocurvature perturbation has no Gaussian part (eq.[4]) and 
thus its cubic term does not vanish unlike the adiabatic case. 
The isocurvature non-Gaussianity therefore can be substan- 
tial even if a ~ 0.067. We neglect the last term B^ss) due 
to the current observational limit on a < 0.067. 



4.2 Bispectra for CMB Temperature Anisotropy 

The angular power spectra for isocurvature perturbations in 
the Quadratic Model are plotted in Fig. [T] for different spec- 
tral indices n a =\, 1.5 and 2. We set a — 0.067 for n„ — 1 as 
in the Linear Model. The power spectrum in the Quadratic 
Model is slightly different from that in the Linear Model due 
to the weak scale dependence of F (eq. [38]). For other spec- 
tral indexes, we set a=0.008 (n a — 1.5) and 0.001 (n CT = 2) 
so that the amplitude of C; around Z=200 becomes roughly 
the same. 

The angular bispectrum from the isocurvature term 
Bsss is approximately given by putting k^ as one of the 
wavenumbers of P aa in the equation (|42[) ; 



,sss 

°hhh 



= 2 



r^r[6f;f (r)&£;7(r)^ ii3 (r) 



+btC(r)b s NLl2 (r)b s L £ s (r) 



J Ll 

+b% Lh (rKf 2 S (r)b s L f/(r)] 



where 
b s,ss 



(r)=~ 



k 2 dkP S s(k)gUk)ji(kr). 



(48) 



(49) 



In order to test the validity of the above approximation, 
we compare the equation (|48|l with full calculation without 
the pole approximation which is used to derive the equation 
(|42p . The angular bispectrum Bsss without the pole ap- 
proxi mation is analytically given as (eq. (C.7) of iKomatsul 
12002] ) 



sss 



p 2 dp 



kUkiP^\(ki, P )g^, 1 (ki)ji[(kir)(-i) h -^ 



%dk2P<-\ (k 2 ,p)g$ l2 (k 2 )j t , for) (-0 



£-1 



kidk 3 g Th (k 3 )j h (k 3 r) + (eye.)], (50) 



where 



I'l'A 



(2l' 1 + l)(2l' 2 + l)(2l + l) 

47T 







h 





X 

and 

P ffCT (|k±p|) 



II 


h 


k 


I'l 


I'l 


1 


h 


1*2 


I 












h I2 h 








1' +i'+i 



im 



(-1) 



i (fc, P )y !m (k)F ; : 



(51) 



(52) 



When the power spectrum is given in power-law form as in 
equation ([7|l. 



k 3 PL\(k, P ) 
2tt 2 



= < 



4:nA a 



4nA 



hi 



ko/ yp / 

when p > k + L^ 1 . 

(53) 

k 



(-l) l Pi~\(k, P ), 



.kj ' 
when p < k — L 1 , 

(54) 



where hi(x) is the expansion coefficient by Legendre function 
Pi(z) as 



hi(y) 



dz{y 2 - 2yz + l) {n "- 4)/2 P,(z). 



(55) 



When n a is equal to 1, hi(y) is analytically given as 2y /(l — 

y 2 )- 

In the left panel of Fig. [3] we compare the bispectrum 
of the adiabatic component with those of the isocurvature 
components with different spectral indexes n CT of 1, 1.5 and 
2. We set / NL = 50 and a = 0.067 for n a = 1, a = 0.008 for 
tict = 1.5 and a = 0.001 for n a — 2 (same as Fig. [TJ. The 
box-size L is set to be 30Gpc. 

Computationally, it is very difficult to evaluate the full 
expression at high I. Thus we check the validity of the pole 
approximations at low £, less than 10. The thick lines repre- 
sent the full calculations of bispectra given in equation (|50p . 
It is found that the isocurvature bispectra approximated as 
the equation (|48|) roughly agree with the full calculations 
within a factor 2 at least for I < 10. The full calculation has 
comparable or larger amplitude for all n a , and a slightly 
steeper slope than the pole approximation. This indicates 
that the isocurvature bispectrum at higher / would become 
larger than the pole approximation and then the proper sig- 
nal may be larger than the pole approximation. The ampli- 
tude of the isocurvature bispectrum is proportional to a 3 ^ 2 
and therefore the effect on the estimation of a is suppressed 
at the power of two-third. 

The isocurvature bispectrum in the Quadratic Model 
depends on the assumed box-size L, which is set to be 30Gpc 
in this analysis. The equation (|45|l indicates that the ratio of 
the isocurvature component in the primordial perturbation 
at k = fco = 0.002Mpc _1 is proportional to a box-size depen- 
dent factor F~ x l 2 (eq. [38]). When L is set to be ten times 
larger (300Gpc), the amplitude decreases by 20% for n a — 1, 
5% for no- = 1-5, and 1% for n a = 2. We find that the over- 
all amplitude of the CMB bispectrum also decreases at the 
same level at / > 10, while additional large-scale power at 
more than 30Gpc slightly increases the amplitude of the bis- 
pectrum at smaller I . The isocurvature bispectrum depends 
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on a 3,/2 and therefore the error in a increases (l/0.8) 2//3 ~ 
16% for jv = 1. 

Using the equation (|30[) . we estimate the isocurvature 
non-Gaussianity in terms of /nl; 



30 



VO.067/ 



3/2 



(56) 



It is found that the isocurvature non-Gaussianity in the 
Quadratic Model can reach /nl ~ 30 given the current 2a 
limit on a. 

The angular bispectrum from the correlation term 
6 Ws = 2 / r 3 *^ (r)6|f (r)6^ ii3 (r) 



+b%LiAr)b C £t(r)b^(r)] 



(57) 



The correlated term b\ \ ' ( is plotted in the right panel of 
Fig. [3] for different spectral indexes n a = l, 1.5 and 2. The 
correlated coefficient cosO^a is set to be 0.1. Using equation 
(|30p . the isocurvature non-Gaussianity with n a — 1 corre- 
sponds to 



/nl = 240 (-^—) 

■' Vo.067/ 



1/2 



COS Ofa 



(58) 



If ^ and a are strongly correlated initially, the correlation 
term can generate substantial non-Gaussianity, while the ef- 
fect of the correlation on Ci is negligible as shown in equation 



5 PERTURB ATI VE FORMULAE OF 
MINKOWSKI FUNCTION ALS 

We adopt the perturb ation formu l ae for Minkowski Func- 
tionals developed by iMatsubaral (|2003h to describe the 
non-Gaussianity of the CMB temperature anisotrop y 
jHikage. Komatsu fc Matsubaral I200II Iffikage et aj]|2008h . 
We separate the analytical formulae of Minkowski Function- 
al into the amplitude and a function of v, which is defined 
as AT/T divided by its standard deviation, as follows; 



Vk{v) = A k Vk{v)- 



(59) 



The amplitude A k is given using the angular power spectrum 
Ci as 



A k = 



">2 



01 



(27r)( fc +D/ 2 U2~ k u k \^/2a 



(60) 
(61) 



where tu k = n h ^ 2 /T(k/2 + 1) gives ujo = 1, Wi = 2, 
UJ2 — t and Wi represents the smoothing kernel deter- 
mined by the pixel and beam window functions and an addi- 
tional smoothing. In our analysis, we add a Gaussian kernel 
Wi = exp [-1(1 + l)<9 2 /2] where S is a smoothing scale. For 
weakly non-Gaussian fields, the function v k (v) can be di- 

( G) 

vided into the Gaussian term v k and the non-Gaussian 
term Avk', 

v k (v) = vi G \u) + Av k (u,f^), (62) 



„( G ) 



Av k (u, / NL ) 



(63) 



-v 2 /2 



±S {0) H k+ 2(v) + ^S w H k (u) 



k ( k ~ !) g(2) 



() 



S w H k -i(v) 



a + O(al) } (64) 



where H n (u) is the n-th Hermite polynomials. The non- 
Gaussian term Av k is characterized by three skewness pa- 
rameters S 1 -^ at lowest order in ao. The skewness parame- 
ters S*-*' are given by the integral of the reduced bispectrum 

as 

3 \ - 2 



-,(oj 



2tto-^ 



E 



(65) 



2<l 1 <l 2 <l 3 



Ml) 



-,(2) 



87T(TQCr 2 



[ii(Ji+l) + Za(ia + l) + la(Js + l)] 



2<1 1 <1 2 <1 3 



xIhhi 3 b W3 W h W h Wi s , 
3 



(66) 



J2 {[h(h + 1) + Uh + 1) - h(h + 1)] 



2<l 1 <l 2 <l'< 



xh(h + 1) + (eye.)} lf ll2l3 b llhh W h W h W l3 , (67) 



where I^^h was previously defined in Eq. (|31|l . 

In Fig. |4j we plot the three skewness parameters 
S 1 ^' from each component of bispectrum in the Quadratic 
Model as a function of 8 3 . Upper panels show the skewness 
of isocurvature components with different spectral indices 
n a =1, 1.5 and 2 in comparison with the adiabatic case. 
The skewness from a mixed component b^ 5 ' is plotted in 
the lower panels. The smaller scale (higher I) information 
in the bispectrum is reflected in skewness parameters with 
a higher number (that is rather than S^) because it 
is weighted towards higher I. We can extract further de- 
tailed scale-dependent information by changing the smooth- 
ing scale 8 3 . 

Fig. [5] illustrates an example of non-Gaussian effects 
on each Minkowski Functional Av k (eq. [BJ) at different 
smoothing scales 8 3 = 10, 20, 40, 70 and 100 arcmin. Here 
we consider isocurvature non-Gaussianity in the Quadratic 
Model with n a = 1 and a = 0.067. It is found that the 
non-Gaussian effect on Minkowski Functionals is 1% or less. 



6 LIMITS ON ISOCURVATURE 

NON-GAUSSIANITY FROM WMAP DATA 

We use the WMAP 5-year data to constrain the non- 
Gaussianity associated with primordial isocurvature pertur- 
bations. We use the linearly co-added maps for Q, V and 
W frequency bands with iV 8 ide = 512. The co-added maps 
are masked with the Kq75 galaxy mas k inclu ding the point- 
source mask provided by I Gold et al.l (|2009h , which leaves 
71.8% of the sky available for the data analysis. The field 
is smoothed with a Gaussian filter at a scale of S . We 
obtain the normalized Minkowski Functional (eq.[64]) for 
the WMAP data using the same procedure described in 
iHikage et all |2008f ). 

We use Bayes's theorem to find a probability of a or 
other combined parameters, together denoted by x, from 
the observed set of Ai/< ota > as follows; 



V(x\Av {ohs) ) oc £{Av {ODB) \x)V(x), 



(obs) | 



(68) 



8 Hikage et al. 



10 

CCC>ru=0.96,a=0 
SSS,n o =r.0,a=0.067 
SSS,n =1 .5,a=0.008 
SSS,n o =2.0,a=0.001 
Full 



100 



1000 




O 



1000 



10" 



10° r (CCS),n a =1"0,a=0.067 
, . (CC5),n„=1 .5,a=0.008 
1CT -(CCS),n a =2.0,a=0.001 



10' - 



£ 10 u 



10" 



-10° 
-10 4 



-10° h 



-10 1 h 

-10 2 h 



10 



100 
■"i — 



f NL =50,|cos9<J=0.1 




1000 



A 



_l I 1 



—I I I ' ' ' ' 1 



1 



10 



100 



—I I 



1000 



Figure 3. CMB angular bispectra of equilateral configurations I 2 (I + l) 2 bm / (2tt) 2 in the Quadratic Model; the isocurvature components 
^flf S (Left) and the mixed components 6>f5 (Right). The power-law index of isocurvature perturbation n CT is set to be 1 (long-dashed), 
1.5 (short-dashed) and 2 (dotted). For comparison, the adiabatic bispcctrum hiS with /nl = 50 is plotted in both panels (thin solid 
lines). The fraction of isocurvature power spectrum a is 0.067 (n CT = 1), 0.008 (n a = 1.5) and 0.001 (n a = 2) defined at ko = 0.002Mpc~ 1 . 
Thick solid lines in Left panels show the full calculations of isocurvature terms (eg. |50J 1 for each n„ at I < 10. The box-size L max is set 
to be 30Gpc. 
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Figure 4. Skcwness parameters S' fe ' (left: k = 0, center: k = 1, right: k = 2) of each component in the Quadratic Model plotted as a 
function of 8 S . The skewness values for a pure isocurvature component are plotted in upper panels, while those for a mixed component 
are plotted in lower panels. The parameters of isocurvature perturbation are n CT =l with a = 0.067 (long-dashed), no-=1.5 with a = 0.008 
(short-dashed), and n tT =2 with a = 0.001 (dotted). For comparison, the adiabatic skewness with = 0.96 and /nl = 50 are plotted 
with solid lines in all panels. 
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arcmin. 



where C is the likelihood function of A?/ obh ' when a non- 
Gaussian parameter has a value x and the probability V(x) 
represents the prior for x. In general, we need to analyze 
all non-Gaussian components together, however we consider 
the cases in which a single non-Gaussian term dominates 
over the other terms for simplicity; in the Linear Model, the 
pure isocurvature term (eq. [A2) ') dominates when there is 
no correlation between <j> and r\. If there exists a weak cor- 
relation, the correlation term of (CSS) (eq. |A8| ) becomes 
important. The other correlation term of (CC<S) (eq. [A7J ) 
dominates when the correlation is strong (cos 6^ ~ 1). In 
these three cases, we set the parameter x to be a 2 /^°', 
a 3,/2 cos ^</)j)/n L ° an d a cos 2 O^-nf^ ^ ■ In the Quadratic 
Model, the pure isocurvature term (eq.[38]) dominates at 
cos #4,0- <C 1. Then we set x to be a (not a 3//2 ). When the 
correlation between <j> and a is strong, the correlation term 
(eq.[57]) is dominant and then x is a 1 ^ 2 cos 2 B$ a . We assume 
a flat prior for all x. We furthermore impose a non-negative 
constraint on the parameters in the Quadratic Model, a and 
a 1//2 cos 2 by definition of a (ea.fTO]). 
The likelihood function is computed by 

-21n£(Au (obs) j:r) oc ^[A^ ob8) - Avf heory) {x)]G^ 

x [A«j ob8) - Avf heoTy) (as)], (69) 

where i and j denote the binning number of threshold values 
v, different kinds of Minkowski Functional k, and smoothing 
scales parameterized with 9„. We choose 18 threshold val- 
ues at an equal spacing in the range of v from —3.6 to 3.6. 
The full covariance matrix dj is estimated from 1000 Gaus- 
sian simulation maps from purely adiabatic perturbations. 
They include the pixel and beam window function, Kq75 
survey mask, and inhomogeneous noise for WMAP 5-year 
maps. Applying our procedure to limit the non-Gaussianity 
from curvature perturbations by putting x as /nl, we obtain 
—63 < /nl < 76 at 95% CL (Hikage et al. in preparation). 

Table[T]lists the mean and one-sigma error of the isocur- 
vature non-Gaussianity in the Linear Model characterized 
by a cos 2 O^f^l , a 3/2 cos O^f^l^, and a 2 f£[°\ Signif- 
icant isocurvature non- Gaussian signals are not found. If 
isocurvature perturbations exist and there is no correlation 
between <f> and r/, the non-linear parameter /^°' is con- 



strained from isocurvature non-Gaussianity for a fixed 01 to 
be 

f^l 0) = (-3300 ± 13000) (a/0.067) ~ 2 . (70) 

If there is a strong correlation between 4> and 77 represented 
by curvaton-type isocurvature perturbations (cos 6^ v = — 1), 
the correlated term (C(S) becomes important and then its 
non-Gaussianity is limited as 

f£l 0) = (4900 ± 43000)(cv/0.0037) _1 . (71) 

Table [5] lists the maximum likelihood value qml, at 
which V(a) has a maximum value, and the 95% confidence 
limit of a when the pure isocurvature term in the Quadratic 
Model dominates. We do not find a significant non-zero value 
of a. The upper limit of a is given by 

a < 0.070 (n a = 1), 
a < 0.042 (n a = 1.5), 

a < 0.0064 (n CT = 2), (72) 

at 95% CL. Our constraint on a for n a = 1 is compara- 
ble to that from the joint limit from WMAP(TT and TE 
spectra)+B AO+SN, which i s sma ller than 0.067 (95% CL), 
obtained bv lKomatsu et all (|2008l 'l. The limits on the corre- 
lated term a 1//2 cos 2 Q<t, a are also listed in Table El Our result 
roughly agrees with iBoubekeur fc Lvthl (|2006h who showed 
the fraction of a Gaussian-squared component of primordial 
curvature perturbation is limited to be less than 0.18, which 
corresponds to a < 0.031, when |/nl| < 100. The difference 
from their analysis is that we calculate a Gaussian-squared 
component of CMB isocurvature bispectra including its full 
radiative transfer function. 



7 IMPLICATIONS FOR AXION 
ISOCURVATURE 

Here we briefly discuss implications of our results in some 
explicit models. Although, for the linear model, the con- 
straint from non-Gaussianity obtained here is not as strong 
as that from the power spectrum, the constraint for the 
quadratic model is severe as we showed in the previous sec- 
tion. Thus we conside r the case of the axion as discussed in 
iKawasaki et al.l (|2008l ). Assuming that Pecci-Quinn symme- 
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try has already been spontaneously broken during inflation, 
the mean value of the axion field can be written as 



fad a 



(73) 



where f a is the axion decay constant and 9 a is the phase 
of the axion. During inflation, the axion field has quantum 
fluctuations 



Sa = 



Hint 

2n ' 



(74) 



where Hi n { is the Hubble parameter during inflation. When 
the average value of the axion is much less the mean-square 
inhomogeneity of a, i.e., when f a a < Hi n t/2n, the density 
fluctuation of a is written as 



Sp a 

pa 



Sa^ 2 
a 



(75) 



Thus power spectrum for isocurvature fluctuations can be 
written as 



k J Pss 
2tt 2 



(n a h 2 



(fi cdm /i 2 ) 2 



(76) 



where Sl a is the energy density of the axion at present. When 
f a 8 < Hi n f/ 2TT, the abundance of the axion is given by 
jTurnedll98a ) 



H ini /27v 
10 12 GeV 



(77) 



Since the total amplitude of primordial perturbations should 
be A to t A cc + A S s = k 3 /2it 2 (P a + P ss ) = 2.4 x 1(T 9 
which is required from WMAP5, by using the constraint on 
a presented in the previous section, a < 0.070, we obtain a 
limit for the Hubble parameter during inflation 



H ini < 1.7 x 10 1 



fa 



10 12 GeV 



GeV, 



(78) 



where we adopt f2 c d m fr 2 = 0.108 which is the mean value for 
a ACDM model from the WMAP5 analysis. The condition 
f a a < H inf /2ir gives 



9 a < 2.7 x 10" 



fa 



10 12 GeV 



(79) 



Although our results obtained in this paper may also 
have implications for other models with isocurvature fluctu- 
ations, a detailed study of this issue will be given elsewhere. 



8 SUMMARY AND CONCLUSIONS 

We have explored the effect of non-Gaussianity from pri- 
mordial isocurvature perturbations on CMB temperature 
anisotropics. Considering the linear and quadratic forms of 
isocurvature perturbations, which are applicable to a wide 
range of theoretical models, we derived theoretical expres- 
sions for bispectra and Minkowski Functionals of CMB tem- 
perature maps with isocurvature non-Gaussianity. We find 
that the amplitude of a quadratic correction f^°^ in the 
Linear Model (a Gaussian variable plus its quadratic cor- 
rection) needs to be of the order of 10 4 to generate CMB 
non-Gaussianity at a level of /nl ~ 10. The isocurvature 
non-Gaussianity in the Quadratic Model (quadratic in a 
Gaussian variable without linear terms) can reach /nl = 30 



while respecting the current upper limit on the isocurva- 
ture contribution to the power spectrum a < 0.067. Isocur- 
vature perturbations provide a possible mechanism to ex- 
plain primordial non-Gaussianity recent ly suggested from 
the observed bispectrum of the CMB (|Yadav fc Wandeltl 
120081 ; iKomatsu etalll2008| ) . 

We give limits on isocurvature non-Gaussianity from 
Minkowski Functionals for the WMAP 5-year data. In the 
Quadratic Model of isocurvature perturbations, we obtain 
a stringent limit a < 0.070 (95% CL) from the non- 
Gaussianity, which is comparable to the current constraints 
from WMAP TT and TE spectra, BAO and SN combined 
a < 0.067 (95% CL). We apply our results to a QCD axion 
isocurvature model and then obtain a limit for the Hubble 
parameter and the phase of the axion. 

We estimate isocurvature non-Gaussianity in the 
Quadratic Model using the pole approximation (eq. [42]). 
The validity can be checked by comparing Minkowski Func- 
tionals with non-Gaussian simulations with a Gaussian- 
squared perturbation. We plan to perform this analysis in 
the near future. 

We employ Minkowski Functionals to characterize non- 
Gaussianity in CMB maps. The application of our work 
to other higher-order statistics is important to utilize non- 
Gaussian information in a more complete manner. Differ- 
ent statistics are sensitive to different aspects of density 
fields and they are affected by possible observational sys- 
tematics (e.g., foregrounds and point sources) in a differ- 
ent way. There is actually some friction between the lim- 
its on /nl from Min kowski Functionals and the bispectrum 
ijHikage et alj 120081 ) ; Minkowski Functionals analysis indi- 
cate a maximum likelihood value around nearly 0, while 
bispectrum analysis favours a more positive /nl around 50 
or more. Complementary analyses with different statistical 
approaches will provide a more robust way to analyze pri- 
mordial non-Gaussianity. 
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APPENDIX A: ADIABATIC, ISOCRUVATURE, 
& MIXED COMPONENTS OF CMB ANGULAR 
BISPECTRUM 

Explicit forms of the angular bispectra from pure 
adiabatic and isocurvat ure mode are calculated as 
(jKomatsu fc SpergelllioOll ) 

b hL = g/NL J r 2 dr[b^(r)b^(r)b% Lh (r) 
+b^(r)b% Ll2 (r)b^(r) 
+&Wr)&$£V)6i'*V)], 



(Al) 



&gf„1, = 2f^> I r J dr[bt^(r)btf 2 "(r)b%u 3 (r) 



(ISO) / 2 , nS, V r, ( r \ h S, VV( ^ h S 





+b S L ^(r)b% Lh (r)b s L ^(r) 






+b% Lh (r)b S Ll T(r)b S Ll T(r)], 


(A2) 


where 










(A3) 


bfr(r) 


= ~ J k 2 dkP vv (k)gUk)Mkr), 


(A4) 


b C NLl (r) 


= \ J ' k2dk 9ii(k)ji{kr), 


(A5) 


b S NLi(r) 


= ~ f k 2 dkgTi(k)ji(kr). 


(A6) 



The angular bispectra from the correlation terms are 
given by 

b (CCS) b ccs b cs< b scc 

hhh ~ hhh hhh hhh 

= 2/^ L SO) J r 2 dr [b^(r)b^(r)b s NLh (r) 

+b^(r)b S NLl2 (r)b^(r) 
+b S NL ll {r)b^{r)b^{r)} 
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Table 1. Mean values and 1 a uncertainties of each isocurvaturc non-Gaussianity parameter in the Linear Model; the pure isocurvaturc 
term c 2 /^°' (eq. lA2l ) and the correlated terms with adiabatic perturbations a 3 / 2 cos ^, ) /^°' (eg, IA8I ) and a cos 2 ^77/^°'' ( e 1- IA7I ). 
The limits are obtained from Minkowski Functionals for WMAP 5-year data at different smoothing scales 8 S and their combination. 



8 S [arcmin] 


a 2 f (lSO) 




a cos 2 


100 


33 ±84 


26 ±62 


91 ± 220 


70 


36 ± 74 


28 ±56 


104 ±210 


10 


71 ± 78 


56 ±61 


256 ± 230 


20 


0±89 


-3 ±66 


73 ± 240 


10 


1±91 


-30 ± 60 


41 ± 230 


10, 20, 40, 70, 100 


-15 ±60 


-18 ± 43 


18 ± 160 



Table 2. Limits on a for the Quadratic Model from Minkowski Functionals for WMAP 5-year data at different smoothing scales 6 S . 
The listed values arc maximum likelihood values «ml an d 95% CL on a. We neglect non-Gaussianity from curvature perturbations and 
their cross-correlation. We impose a non- negative condition on a following from its definition in equation 11101 . 



0s [arcmin] 


n a 


= 1 


n a 


= 1.5 




= 2 


"ML 


95% CL 


«ML 


95% CL 


«ML 


95% CL 


100 


0.037 


< 0.115 


0.022 


< 0.100 


0.005 


< 0.047 


70 


0.036 


< 0.105 


0.029 


< 0.079 


0.013 


< 0.030 


10 


0.051 


< 0.112 


0.010 


< 0.057 





< 0.014 


20 





< 0.101 





< 0.048 





< 0.0085 


10 





< 0.093 





< 0.047 





< 0.0071 


10, 20, 40, 70, 100 





< 0.070 





< 0.042 





< 0.0064 



2 , observational constraints and thus they are neglected in the 

+-jnl r dr x 



5 J following analysis. 

[bi$+ir) {b^(r)b< NLls (r) + b^(r)b< NLh (r)} 

(r) {bfC'(r)b% Ll2 (r) + b s L ^(r)b% Lh (r)}] , 

(A7) 



b ((SS) _ b css , b scs , b ssc 

= ^/nl J r 2 dr[b S ^{r)btf 2 -{r)b% Ll3 {r) 

+b S L it V ( r ) b NU 2 (r)b S L ^(r) 
+b% Ul (r)b S L ^(r)b S L ^(r)] 

(r) {bt^(r)b% M3 (r)+b s Ll ^(r)b s NLl2 (r)} 
+b ( L f 2 V (r) {bf^(r)b S NLh (r) + b S L ^(r)b S NLh (r)} 
+b ( L ? 3 v (r) {b S Ll ^(r)b s NLl2 (r) + b s L ™(r)b s NLh (r)}] , 

(A8) 

where 

b C £f>(r) = \ f ^dkP^g^kMkr), (A9) 
>fi*"(r) = 1 j k 2 dkP H (k)gUk)ji(kr). (A10) 



The terms proportional to /nl in bj^f'j' an d are very 

small relative to the adiabatic bispectrum under the current 
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Table 3. Same as Tablc[2]but for the parameter a 1 / 2 cos 2 where the correlation term hi the Quadratic Model dominates. We 

also impose a non-negative condition on q 1 / 2 cos 2 6^, a . 



(a 1 / 2 cos 2 e 0CT ) ML 95% CL (a 1 / 2 cos 2 0ct ) M l 95% CL (a 1 / 2 cos 2 <V)ml 95% CL 



100 


0.018 


< 0.095 


0.018 


< 0.12 


0.010 


< 0.089 


70 


0.021 


< 0.092 


0.021 


< 0.11 


0.014 


< 0.073 


■10 


0.045 


< 0.12 


0.040 


< 0.14 


0.017 


< 0.075 


20 


0.012 


< 0.098 


0.026 


< 0.11 


0.013 


< 0.041 


10 





< 0.071 


0.003 


< 0.063 


0.0026 


< 0.019 


10, 20, 40, 70, 100 





< 0.049 





< 0.039 





< 0.012 



